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1 Introduction 

Various extensions of the Standard Model (SM), such as left-right symmetric models and 
supersymmetric models with R-parity violation, predict non-standard interactions (NSIs) 
of the neutrinos with other fermions [1-11]- The NSIs in those models are generated via 
the exchange of new massive particles and at low-energies can be described by effective 
four fermion interactions of the form 

£nsi = -2V2 GFe{/^^{i7^j^'PLUf,) {7 j^Pc f) , (1.1) 

where Gp is the Fermi constant, a and /3 are neutrino flavor indices, / and /' label light 
SM fermions, C = L or i? is the chirality of the projection operator Pq, where -P^/ij = 
(1=F 75)72, and the dimensionless number e^^ parametrizes the strength of the interaction. 

In a previous publication from 2002 [12], Raghavan, together with Berezhiani and 
Rossi, discussed the potential of the Borexino detector in placing constraints on the flavor- 
diagonal NSI parameters 

eaL = £aa > ^aR = Saa > a = 6 Or T , (1.2) 

via the measurement of the electron recoil spectrum in fQ,e elastic scattering. There, it 
was argued that due to the mono-energetic nature of the ^Be solar neutrinos, Borexino 



- 1 - 



would be able to place stronger constraints on £eR and e-tR than would be possible at 
Super-Kamiokande (SK) and the Sudbury Neutrino Observatory (SNO) where the observed 
neutrinos are the neutrinos with a continuous energy spectrum. The a = /i case was 
not considered since the couplings and e^/j were already constrained to the level of 
k/xL/ijI < 0.03 (at 90% C.L.) [13] by the ly^e scattering experiment CHARM II [14]. 

Today, a decade later, with the Borexino experiment running smoothly and having 
accumulated more than 153.6 ton-years of data, it is now possible to actually extract the 
constraints on the flavor-diagonal NSI parameters discussed by Raghavan et al. Constraints 
on the same parameters from various other experiments are also available for comparison, 
some of which are quite new. Bounds from the solar neutrino experiments (SK, SNO, 
etc.) and KamLAND can be found in Refs. [13], [15], and [16]. Bounds from reactor 
and accelerator experiments have been compiled in Ref. [17], which includes the bounds 
from e^e~ — >• vv'^ measured at LEP [18], i^gC scattering measured at LSND [19], and Vge 
scattering measured at the reactor experiments Irvine [20], Rovno [21], and MUNU [22]. 
See also Refs. [23, 24]. Bounds have also been placed using atmospheric neutrinos [25] and 
MINOS [26]. New bounds from the TEXONO reactor neutrino experiment can be found 
in Ref. [27]. 

In this paper, we place constraints on the NSI parameters £eL/R a-^d £tL/r using 
the Borexino 153.6 ton-year data. We include in our analysis the uncertainties in the ^Be 
solar neutrino flux, and backgrounds from ^^Kr and ^^^Bi /3-decay. Taking these systematic 
uncertainties and backgrounds into account are particularly important since they can mimic 
non-zero values of the parameters in question. Indeed, as will be discussed in detail in the 
following, we find that unless the solar-neutrino flux- uncertainty and the ^^Kr /3-decay 
background are reduced, increased statistics will not improve the bounds beyond what can 
be extracted from current Borexino data. The parameters £^l/r are not considered since 
the CHARM II bound stiU stands. 

This paper is organized as follows. In section 2, we discuss how the electron recoil 
spectrum in f^e scattering can be used to constrain the parameters £aL/R^ (a = e or r) and 
how the above mentioned systematic uncertainties and /3-decay backgrounds can interfere 
with the extraction. Section 3 discusses how well systematic uncertainties and backgrounds 
are understood in the Borexino experiment, section 4 discusses the uncertainty in the ^Be 
neutrino flux in the Standard Solar Model, and section 5 discusses the uncertainties in the 
neutrino oscillation probabilities. In Section 6 we look at the dependence of the electron 
recoil spectrum on the NSI parameters and backgrounds in detail, and how one can mimic 
the other. Results of our flt are presented in section 7, and section 8 concludes. 



2 Neutrino-Electron Elastic Scattering 

In the SM, the interaction between neutrino flavor a {a = e, r) and the electron is 
described at low energies by the effective four fermion interaction 



SM 



GpiUa-f'^PLl^a) galinfiPie) + gaRie'y^^PRe) 



(2.1) 
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The coupling constants at tree level are given by gaR = sin^ 9w and gaL = sin^ 6w ± | , 
where the lower sign applies for a = fj, and r (from Z exchange only) and the upper sign 
applies for a = e (from both Z and W exchange). Radiative corrections are small and can 
be ignored.^ 

The presence of flavor-diagonal NSIs, SaL/R^ "^ill shift the coupling constants in the 
above expression to 

9aL QaL = QaL + ^aL , QaR gaR = gaR + £aR ■ (2.2) 

This interaction between electron and neutrino, with possible shifts in the coupling con- 
stants, can be observed via the elastic scattering of a neutrino of flavor a off of an electron 
at rest, which has the differential cross section [30] 



dT TT 



gaL + gaR \ ^ " ^ ) " gaigaR^^ 



(2.3) 



Here, rrig is the electron mass, Ey^ is the initial neutrino energy, and T is the kinetic energy 
of the recoil electron which has the range 

E 

< T < Tmax = • (2.4) 

If the incoming neutrino beam is mono-energetic, no convolution of Eq. (2.3) with the 
neutrino energy spectrum is necessary. 

The "^Be solar neutrinos are produced via the K-shell electron capture processes [31] 

^Be + e- ^ r^^ + ^^ (2.5) 
[^Li*(0.48) + z.e (10.4%), 

yielding mono-energetic neutrinos of energy 0.862 MeV and 0.384 MeV, respectively. Borex- 
ino is sensitive to the 0.862 MeV component, for which the maximum recoil energy is 
Tlnax = 0.665 MeV. This component arrives at the Earth as a superposition of the three 
neutrino flavors due to the MSW effect and vacuum oscillation. If we denote the sur- 
vival probability of Ve in this component at the Borexino detector as Pee, then the total 
neutrino-electron scattering cross section there will be the combination 



day{T) 
dT 



dT ^ ""'V dT dT 



(2.6) 



^In the notation of Ref. [28], the coupling constants for the a = n case, including radiative corrections, 
are expressed as 

gi^L — ^ — Pi^e + H^eSz j , 9^/? — ^ — P''''{^2 

where s% is the MS value of sin^ 9w, and p^e and K^e denote process specific corrections. At zero-momentum 
transfer, these are p^/e = 1.0128 and k^e = 0.9963 [28]. Since the deviations of these parameters from one 
are small compared to the sensitivity of Borexino, they, and any flavor dependence that may exist for the 
a = e, T cases, can be ignored. In our analysis, we use the MS value of (0.862 MeV) = 0.2386 [29] for 
sin'^ 6w, and the tree-level expressions for the couplings. 
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where C23 = cos^ ^23 and S23 = sin^023- Thus, the measurement of the T-dependence of 
the ^Be solar neutrino elastic scattering event rate will let us constrain the values of gai 
and gaR for all flavors. 

An actual detector, however, cannot measure the recoil electron energy to arbitrary 
precision, and one must take its finite energy resolution into account. If we denote the 
probability of detecting energy T4 for an electron with kinetic energy T by RiT^^T), the 
differential cross section as a function of the actual detection energy Ta is given by 

(2.7) 



dTA 



R{Ta,T) "^^^-i^^ dT 



dT 



For Borexino, we take the energy resolution function R{Ta,T) to be a gaussian with a 
T-dependent standard deviation 



R{Ta,T) 
where cr(T) is given by [32] 

a{T) 



1 



27ra(r) 



exp 



MeV; 



1/2 



■2[a(r)P 



ao = 50keV. 



(2i 



(2.9) 



This will 'blur-out' the shape of the energy spectrum somewhat, smoothing out the Compton- 
like edge of da^{T)/dT, but still retain its basic overall shape. 

Thus, for an incoming 0.862 MeV ^Be neutrino flux of ^t'^^^ and number of electrons 
Ne in the fiducial volume of the detector, the number of recoil electrons detected with 
energy in the bin Ti < Ta < T2 per unit time is given by 

dN{Ti,T2) 



dt 



Ti 



A 



dTA 



Ti 



dTA 

dTA 
^2 d a.^TA) 
dTA 



dTA 
dTA 



(2.10) 



By measuring this spectrum, Borexino can constrain both eeL/R ^tL/r- However, the 
precision of those constrains will depend on the uncertainties in the prefactors [Nf,^^'^^'^ Pf>e] 
and [iVe$7gg^(l — Pee)s2^], which are still quite significant. 

The measurement of the recoil electron energy spectrum in Borexino is further compli- 
cated by the fact that it is impossible to distinguish between electrons from fQ,e scattering 
and those from /3-decay of radioactive nuclei. The most significant /? backgrounds in Borex- 
ino are those from the decays 

^^Kr — > s^Rb + e- +Ve (Q = 0.687 MeV, i^/^ ^ 755 yg^rs, 99.57%) , 

2iOBi — > '^^^Fo + e- +Ve (Q = 1.161 MeV, = 5.012 days, 100%) . ^' ' 

As will be shown later, the electrons from ^^Kr decay are particularly problematic. Of the 
decay products, ^^Rb is stable while ^^"^Po subsequently undergoes a-decay 



210 



Po 



206pb + ^He 



(Q = 5.3 MeV, t^/^ = 138 days, 100%) 



(2.12) 
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This Q-particle can mimic /3-particles in the relevant energy range, but this particular 
background can be removed reliably using pulse shape analysis [33]. ^'^^Pb is stable. 

In the next section, we will review the properties of the Borexino detector, and how 
well these /3-decay backgrounds are understood. 

3 The Borexino Detector and beta-decay Backgrounds 

Borexino is a real-time solar neutrino detector designed to measure the 0.862 MeV mono- 
energetic ^Be solar neutrinos. It is situated in Hall C of the Laboratori Nazionali del Gran 
Sasso (LNGS), Italy, below 1400 meters of rock (3800 meters water equivalent) where the 
muon flux is suppressed by a factor of ~ 10^ compared to the Earth's surface. Borexino's 
spherical vessel is filled with ~278 tons of liquid scintillator (pseudocumene doped with 
2,5-diphenyloxazole as a wavelength shifter), with the fiducial volume consisting of the 
central 100 tons. In the period from May 16, 2007 to May 8, 2010, Borexino had 740.7 live 
days of data taking, corresponding to 153.6 ton-years of fiducial exposure. [34-36] 

The neutrinos are detected via their elastic scattering off of electrons in the scintillator. 
The scintillation light from the recoil electrons spread isotropically from the event location, 
and are detected by an array of ~2200 photomultiplier tubes (PMTs) mounted onto a 
stainless steal sphere, looking into the detector volume. A neutrino event is identified by 
multiple PMT hits within a Trigger Time Window (TTW) of 60 nanoseconds, with the 
total number of photoelectrons collected exceeding a given threshold. This triggers the 
recording of all PMT hits for 16 microseconds, from which the electron recoil energy Ta 
is reconstructed. The location of the event within the detector is determined from the 
difference in arrival times of the TTW photons to different parts of the PMT array. For a 
Ta ~ 1 MeV event, this system can identify its location to a precision of about 10 cm [37]. 

The chemical composition of pseudocumene is C9H12 with atomic mass 120.19, and 
66 electrons per molecule. Thus, the number of electrons within the fiducial volume of 
100 tons of pseudocumene can be calculated to be iVg = 3.307 x 10^^. However, due to the 
limited event-position resolution of Borexino discussed above, there exists an uncertainty 
in the fiducial volume of % which propagates directly into an uncertainty in A'^e [36] . 

As discussed in the previous section, the electrons from ®^Kr and ^^'^Bi /3-decay present 
significant backgrounds to the measurement of the ^Be signal. Due to the Q-values of these 
decays being close to the ^Be solar neutrino energy, their /3-decay spectra overlap with the 
^Be recoil electron spectrum as shown in figure 1. Of the other backgrounds present, that 
due to ^^C, which /^-decays with Q- value of 0.1565 MeV, is large but occupies a much lower 
energy range [38]. The background due to cosmogenic ^^C [39], which /^"'"-decays with 
Q- value of 0.96 MeV, is at larger energies. A fit to the measured count rate assuming SM 
interactions yields (TABLE I of Ref. [36]) 



^Be 

85Kr 

210Bi 



46.0 ± 1.5(stat)li;g(syst) counts/(day • 100 tons) 
31.2 ± 1.7(stat) ± 4.7(syst) counts/(day • 100 tons) 
41.0 ± 1.5(stat) ± 2.3(syst) counts/(day • 100 tons) 



(3.1) 
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showing that the count rates of the ^^Kr and ^^'^Bi backgrounds are of the same order 
of magnitude as that of the ^Be signal. Thus, the distortion of the shape of the count 
rate spectrum due to these backgrounds are significant. Furthermore, the shapes of the 
®^Kr and -^^"Bi /3-spectra are such that the said distortions are similar to those due to 
non-zero values of the NSI parameters. Thus, in order to extract the bounds on the NSI 
parameters from the count rate data, it is crucial that we can determine the ^^Kr and ^^"Bi 
backgrounds independently. 

For ^^Kr, this is possible by utilizing the following decay chain which constitutes 0.43% 
of ®^Kr-decay : 

85Kr »5^Rh + + z7e (Q = 0.173 MeV) , 

^5™Rb s^Rb + 7 (Q = 0.514 MeV, t^/^ = iq^S) . ^ ' ' 

Delayed coincidence measurements of the (3 and 7 from this decay-chain has yielded the 
total count rate of [40] 

^^Kr : 30.4 ± 5.3(stat) ± 1.3(syst) counts/(day • 100 tons) , (3.3) 

which is consistent with the fit value listed in Eq. (3.1). We will use this direct measurement 
value constraint in our analysis. 

^^"^Bi is a pure /3-emitter produced at the end of the ^^^Rn decay chain: 

222^^ ^ 218p^ ^ 214pb ^ 2143. ^ 214p^ ^ 210pb ^ 2103^ ^ 210p„ ^ 206pb _ 

(3.4) 

It was proposed in Ref. [41] that the quantity of ^^"^Bi in liquid scintillation detectors may 
be determined directly from the quantity of its decay product ^^"Po, which in turn can 
be determined by measuring its a-decay rate. Unfortunately, this method cannot be used 
due to temporal instabilities of the Borexino data. However, in the narrow energy range 
above the ^Be shoulder and below the lower end of the ^^C /3~^-decay spectrum, aka the 
"^Be valley," the count rate is dominated by ^^'^Bi decay. The count rates from the other 
two components that contribute in this valley, namely the CNO and pep solar neutrinos, 
are too low compared to that of ^^'^Bi to be measured with statistical significance. Indeed, 
the contributions of CNO and pep neutrinos are small throughout the energy range of our 
analysis, with count rates comparable to the statistical errors of the /3-decay backgrounds 
[42]. Therefore, though the precision is limited by statistics, fitting to the count rate in the 
valley can constrain the ^^"Bi background independently of the fit to the NSI parameters. 



4 Uncertainty in the Standard Solar Model Neutrino Flux Prediction 

The precisions of the neutrino fiuxes predicted by the Standard Solar Model (SSM) depend 
on the precisions of the physical input parameters. In the past, the "^Be neutrino fiux 
prediction suffered from a large uncertainty (12% in Ref. [43] from 2004) stemming mostly 
from a large uncertainty in the cross section of the reaction 

^He + ^He ^ ^Be + 7 (4.1) 
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at energies relevant to reactions in the Solar core. This cross section had been measured 
by shooting an a-beam into a gas ^He target, and counting the number of produced ^Be 
by either detecting the 7s in the above reaction (prompt method), or by measuring the 
amount of accumulated '^Be later by counting the 7s from the decay of '^Li*(0.48) produced 
in the second reaction of Eq. (2.5) (activation method). The large uncertainty was due 
to a disagreement between the results obtained via the prompt and activation methods. 
Extrapolating the results down to the relevant energies lead to an additional uncertainty. In 
the last decade, however, these uncertainties have been reduced considerably by the LUNA 
experiment at the LNGS by a resolution of the discrepancy between the two methods, and 
the measurement of the cross section at lower energies than before [44, 45]. As a result, 
the error in the ''Be neutrino flux due to this particular uncertainty has been reduced from 
8% to 2.8% [46]. 

Unfortunately, solar metalicity has come into the picture as a new source of large uncer- 
tainty. The abundance of various elements inside the Sun is inferred from the measurements 
of their abundances in planets, chondrites^, and other bodies in the Solar system, and from 
the analysis of absorption lines in the Solar photosphere spectrum. The abundances in- 
ferred from the latter depend on the hydrodynamic model used for the photosphere. A 
determination based on a ID model of the solar atmosphere published in 1998 [47] (GS98) 
has been used as input to the SSM for many years. A new determination based on a 
state-of-the-art 3D model, and other improvements, was published in 2009 [48] (AGSS09). 
The inferred values of the metalicity at the solar surface is (Z/X)q = 0.0229 in GS98, but 
(Z/X)0 = 0.0178 in AGSS09.3 The SSM predictions of the total ^Be neutrino flux (sum 
of the 0.862 MeV and 0.384 MeV components) for the two cases are [49] 

= 5.00 (1 ± 0.07) X 10^ cm-^s'i (GS98) , 

= 4.56 (1 ± 0.07) X 10^ cm-^s-i (AGSS09) . (4.2) 

The 7% uncertainty in both cases are due to the uncertainties in the cross sections of the 
reactions ^He + ^He ^He + 2p (2.5%) and ^He + ^He ^ '^Be + 7 (2.8%), opacity (3.2%), 
diffusion (2%), quoted uncertainties in the abundances (2%), and various other sources [50]. 
As can be seen, the central value of the flux prediction can differ considerably depending 
on which set of abundances are adopted. However, various predictions of the SSM based 
on the AGSS09 data set disagree with solar properties inferred from helio-seismology, while 
those based on GS98 show excellent agreement. (See e.g. Ref. [51] and references therein.) 
Thus, we will simply use the GS98 based flux for our analysis, which was also the flux used 
by Borexino to constrain Pee [36]. The 0.862 MeV neutrinos constitute 89.56 it 0.04% [31] 
of the total flux, so the GS98 prediction of ^7^^'^ wih be 

$0j862 ^ 4 43 ^ Q_Q7^ ^ -^q9 cm-^g-i (GS98) . (4.3) 
kind of meteorite. 

^X, Y, and Z respectively denote the mass fractions of hydrogen, heUum, and metals (elements other 
then H or He). 
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l(j bound, from R.ef. [53] 


l(j bound from R,ef. [54] 


Am|i (lO-^eV^) 


7.62 ±0.19 


7 r.A+0.26 
'■^^-0.22 


A/vTo2 /in— 3^A/"2\ /'ATXJ^ 
l\m^-^ eV ) ^iNrlj 

13 V / V / 


9 cq+O.OS 

^•"-•"J-o.io 

9 4^+0.10 


9 /iq+0.06 
^•^-S-O.IO 

2 42+0-07 


sin^6'i3 (NH) 
sin^^ia (IH) 


027+0-003 


0241 ± 0025 
0.0244^0.0^1 


sin^ 023 (NH) 
sin^ 023 (IH) 


n 40+0.08 

'-'•^^-0.05 


0.386lH^t 

Q 009+0.039 

u.oyz_Q 022 


sin^ 012 


320+0-015 

o.O^U_o 017 


307+0-018 


p 


n 5fi^+o-oi3 


t^74+0-0i2 

'^•^'*-0.013 



Table 1. Global fit values of the neutrino mixing angles and mass-squared differences from Refs. [53] 
and [54], and the corresponding ly^ survival probability Pee for the ^Be neutrinos. NH and IH stand 
for normal and inverted hierarchies. The agreement is good for all parameters except sin^023, the 
value of -which differ considerably for both hierarchies. 

5 Uncertainties in the Mixing Angles 

As discussed in section 2, -we need to understand the uncertainties in the factors [Ne^r'^^'^ Pee] 
and [A''e$7gg^(l — Pee)s2s\ ill Eq. (2.10) to constrain s^l/r (" = e, r). In the previous sec- 
tions, -we have seen that A^e has an uncertainty of %, -while ^r^g^ has an uncertainty 
of ±7%. Let us now look at the uncertainties in the remaining factors. 

The fe survival probability Pee is dependent on the neutrino energy Ei,. At Ei, = 
0.862 MeV, matter effects are negligible and -we can approximate (see, for instance, Ref. [52]) 

Pee ~ 1- ^Sin2(20i2) . (5.1) 

The best-fit values of the neutrino mixing angles from t-wo recent global fits, Refs. [53] and 
[54], are listed in Table 1. These analyses are more recent that -what -was used in the 2012 
Revie-w of Particle Properties [55] . The agreement in the values of sin^ 0i2 from the t-wo 
fits is good, and rounding up the resulting Pee to the second decimal place yields 

Pee = 0.57 ±0.01 (5.2) 

for both. We -will use this value for Pee in our analysis.^ Thus, the uncertainty in the factor 
[A^e^7Be^-^ee] is dominated by that in ^I^rgf^, and we assign to it a ±7% uncertainty. 

The value of sin^ 023 is problematic. While the value of Ref. [53] is consistent -with 
the maximal mixing, sin^ 023 = 0.5, Ref. [54] claims a deviation. The claimed precision of 
the values are also quite different. Since it is beyond the scope of this paper to settle this 
controversy, we fix sin^ 023 to the maximal mixing value of 0.5 for the purpose of this paper, 
and ignore the uncertainty in the factor (1 — Pee)s23- Thus, the factor [Ae<I'7Be(l — -Pee)'S23] 
is assigned a ±7% uncertainty also. This simplification should be kept in mind as a caveat 
in interpreting the robustness of our bounds on e^j^j^. 

*The value based on the Borexino data assuming Eq. (4.3) is Pee ~ 0.51 ± 0.07 [36]. 
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Figure 1. Expected Borexino event spectrum for a fiducial exposure of 153.6 ton-year. The 
numbers shown are for lOkeV wide bins, with associated statistical errors. The left-hand panel 
depicts the total number of events (''Be + ^^Kr + ^^^Bi) for different choices of EeL and EcR as shown 
in the legends. The contributions from ^^Kr and ^^"Bi, which are unaffected by the NSI parameters, 
are also shown separately. The right-hand panel portrays the same for different choices of and 
Et-r- Note that we have taken the sample non-zero values of the NSI parameters to be fairly large 
to exaggerate their effects. 



6 Dependence of the Event Spectrum on the NSI Parameters 

Before we proceed to our analysis of the Borexino 153.6 ton-year data, let us take a look 
at how the expected event spectrum depends on the NSI parameters to obtain a feel on 
which features of the spectrum are relevant in constraining which parameter. 

In Fig. 1, we show the expected 153.6 ton-year Borexino event spectrum in the energy 
range 0.29 MeV < Ta < 0.80 MeV, in lOkeV wide bins, for several choices of the NSI 
parameters. This energy range is dominated by the ^Be, ^^Kr, and ^^''Bi events, and 
all other solar neutrinos and background contributions can be neglected. The ^Be signal 
has been calculated with inputs sin^ 6]y = 0.2386, Pee = 0.57, s^s = ^^23 — O-^; ^^'^ 
$7g^^ = 4.48 X 10^cm~^s~^, as discussed in previous sections. The total number of ^Be 
counts in the shown range is 14350 for the SM case (all NSI parameters set to zero). The 
®^Kr background has been fixed to the rate in Eq. (3.3), and the ^^*^Bi background to the 
rate in Eq. (3.1). The total number of counts in the range is 5813 (10057) for the ^^Kr 
(^^*^Bi) background. 

It is quite evident from Fig. 1 that the left-handed couplings, eaL(a = e,r) affect the 
overall normalization, whereas the right-handed couplings, £aR{oi = e, r), cause changes in 
both shape and normalization. Thus any uncertainty in the normalization of the ''Be signal 
can severely deteriorate the sensitivity to £aL , and also mimic the presence of SaR to some 
extent. The ®^Kr background starts around the Compton-like edge of the ^Be signal, and 
its presence affects the slope of the total event spectrum below this edge in a way similar 
to non-zero values of SaR- Thus, the uncertainty in this background can be expected to 
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-0.2 -0.1 0.1 0.2 -0.1 0.1 



^eL I ^eR ^eL ^ ^eR 

Figure 2. Correlations between the NSI parameters and the normalizations of the ^Be signal and 
^^Kr background. The left-hand panel shows the 95% C.L. contours (2 d.o.f., Ax^ — 5.99) in the 
AiVBe-eeL,ij plane, while the right-hand panel shows the same in the ANKi-£eL,R plane. The area 
inside the contours are allowed. All implicit parameters are set to their reference values used in 
Fig. 1. 



deteriorate the sensitivity to SaR- Note also that the events above the ^Be Compton-Iike 
edge, 0.7 MeV < Ta, are dominated by the ^^'^Bi background as discussed previously, so 
this background will be constrained by the events in this region when performing a fit. 

To see how the uncertainties in the normalization of the ^Be signal and the ^^Kr 
background could affect the bounds on the NSI parameters, we perform the following 
analysis. Let A^?^ be the expected number of events in the i-th energy bin. Construct the 
by 

X (sacANx) = ^ N'^{0,0) ' « = e or r , (6.1) 

where C = L or R, and X = Be or Kr. AA'^bc and AA^Kr are respectively the percentage 
changes in the normalizations of the ^Be signal and '^^Kr background from their reference 
values. Implicit parameters are all fixed to the reference values we used above. This 
quantifies the ability of Borexino to distinguish between the {SaCi ^-^x) = (0, 0) and non- 
zero cases. Using this for the case a = e, we plot the 95% C.L. contours (Ax^ = 5.99) 
in the ANx-£eC planes in Fig. 2. 

The left-hand panel of Fig. 2 shows the correlations between £eL/R and AA'^Be- For 
EeL, we can see that it has a strong negative correlation to AA^bc, as was expected from 
our discussion above. Thus, a large uncertainty in AA^bc would lead to a large uncertainty 
in EeL- ^eR, OH the other hand, is only weakly correlated with AA^Be, again as expected. 

The right-hand panel of Fig. 2 shows the correlations between £eL/R and AA'xr- For 
EeL, we see that it is only weekly correlated with AA^'xr as expected. EeR, however, is 
also only weakly correlated with AA^xr somewhat contrary to expection. Consequently, 
a reduction in the uncertainty of AA'^Kr will not lead to any significant reduction in the 
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uncertainty of EgR. In the following, we will find that the reduction of the Kr background 
itself, and not just its uncertainty, is necessary to improve the bound on EeR- 

7 Analysis Results 
7.1 Method of Analysis 

Let us now proceed to our analysis of the Borexino 153.6 ton-year data. As in the previous 
section, we consider events in the energy range 0.29 MeV < Ta < 0.80 MeV, which we 
divide into 10 keV bins. 

Let the number of measured counts in the i-th bin be N^^"^, and its theoretical value 
A^*'^(A), where A denotes the fit parameters that will be varied: 



AA^Bei AA^Kr and AAsi respectively denote the percentage change in the ^Be, ^^Kr, and 
^^"^Bi event normalizations from their reference values. The first is varied to account for 
the common ±7% uncertainty we assigned to the three prefactors in Eq. (2.10), while the 
second for the ±18% uncertainty in the *^Kr background, c.f. Eq. (3.3). The is then 
defined as 



No prior constraint is imposed on AN^i, which will be left for the fit to determine. 

Since we do not have access to the raw 153.6 ton-year Borexino data, we reconstruct 
the counts N^^"^ from the fits listed in Eq. (3.1) as a sum of ^Be, '^^Kr, and -^^"Bi events. 
Contributions from other solar neutrinos and backgrounds are neglected. The theoretical 



Using this Ax^, we first place constrains on the four NSI parameters one at a time, keep- 
ing the other three NSI parameters zero, while marginalizing^ over the three normalization 
parameters. Then, we place constraints on pairs of NSI parameters, one flavor at a time, 
keeping the other flavor pair zero, while again marginalizing over the normalization pa- 
rameters. 

7.2 One NSI parameter at a time limits 

The dependence of Ax^ on one NSI parameter, with the other three fixed to zero, and 
after marginalization of the three normalization parameters, is show in Fig. 3 for all four 
choices of the NSI parameter. The corresponding 90% C.L. limits (1 d.o.f, Ax^ = 2.71) are 
listed in Table 2, together with the global fit values from Ref. [17] for comparison. We can 

^For a discussion on the 'marginalization' procedure see, for instance, Refs. [56] and [57]. 



A = {SeL, SeR, £tL, ^tR, AAbg, AiVKr, AiVBi} • 



(7.1) 




(7.2) 
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One parameter at-a-time One parameter at-a-time 




-0.2 -0.1 0.1 0.2 -1 -0.5 0.5 1 

^eL I ^eR ^xL ' ^xR 



Figure 3. Left-hand panel depicts the one NSI parameter at a time hmit for e^L and e^r- 
Right-hand panel shows the same for Etl and SrR- The corresponding limits at 90% C.L. (1 d.o.f, 
Ax^ = 2.71) are listed in Table 2. 













This work 


[-0.046, 0.053] 


[-0.206, 0.157] 


[-0.23, 0.86] 


[-0.968, 0.716] 


Global limits [17] 


[-0.03, 0.08] 


[0.004, 0.151] 


[-0.5, 0.2] 


[-0.3, 0.4] 



Table 2. The 90% C.L. limits on the flavor-diagonal NSI parameters CeL, Eefi, £tL and £rR 
based on 153.6 ton-years of Borexino data. In each case, only one NSI parameter and the three 
normalization parameters are allowed to float, while the remaining three NSI parameters are fixed 
to zero. The three normalization parameters are marginalized to obtain these bounds. The second 
row lists the global bounds from Ref. [17] for comparison. 



see that the bounds based on Borexino data alone is already competitive with the global 
fit to reactor -|-acceIerator data. 

Since July 2010, the Borexino experiment have undertaken a series of purification 
campaigns to reduce the radioactive backgrounds. The method of Nitrogen stripping has 
been quite successful in removing the ^^Kr background by roughly 70% [58]. In 2012, with 
these lower backgrounds, the Borexino experiment has entered into its Phase II run and 
continues to accumulate more data. 

Let us quantify the expected improvements on the NSI bounds due to improved statis- 
tics and the reduction of the ^^Kr background. In the leftmost panel of Fig. 4, we show 
how the bounds on the electron NSI parameters SgL and EeR will be affected by an increase 
in the total fiducial exposure while keeping all other assumptions the same. The vertical 
dot-dashed line indicates the current fiducial exposure of 153.6 ton-years. The leftmost 
panel of Fig. 5 shows the same for the tau NSI parameters £rL and EtR- We can see from 
these panels that increased statistics will not improve the limits on the left-handed coupling 
at all, while the improvements in the right-handed couplings are quite modest, indicating 
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Background [%] 
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Figure 4. One NSI parameter at a time limits on and SeR at 90% C.L. (1 d.o.f) as a function 
of total fiducial exposure (left panel), total amount of ^^Kr background (middle panel), and la- 
uncertainty on ^Be signal rate (right panel). 



that the current uncertainties in the NSI parameters are mostly due to backgrounds and 
systematic uncertainties. 

In the center panel of Fig. 4, we show the impact of a reduction in the ^^Kr background 
on the bounds on EeL and SeR, with a fixed fiducial exposure of 153.6 ton-years. The 
reduction is expressed in percentages compared to Eq. (3.3). The center panel of Fig. 5 
shows the same for e^-l and ErR- As discussed in section 6, we expect the bounds on 
the left-handed couplings, which change the ^Be signal normalization, to be little affected 
since the ^^Kr background mostly changes the slope of the spectrum, and indeed the figures 
confirm this expectation. The bounds on the right-handed couplings, on the other hand, 
can be tightened. If the ^^Kr background is reduced down to 10%, the bounds on eg/j will 
shrink by a factor of ~2, while those on ErR will strink by a factor of ~ 1.5. 

In the rightmost panels of Figs. 4 and 5, we show the impact of a reduction in the 
uncertainty of the ''Be signal normalization on the bounds on EeL and EeR, and ErL and 
ErR, respectively. All other assumptions, including the fiducial exposure of 153.6 ton-years, 
are kept the same. The vertical dot-dashed lines show the current la uncertainty of ±7%. 



1 1 1 1 1 1 1 


' ' ' ' 1 1 1 
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Background [%] 



90% C.L. (1 d.o.O 
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Figure 5. One NSI parameter at a time limits on ErL and ErR at 90% C.L. (1 d.o.f) as a function 
of total fiducial exposure (left panel), total amount of ^^Kr background (middle panel), and la- 
uncertainty on ^Be signal rate (right panel). 



- 13 - 




Figure 6. Allowed regions in the EeL-Eefl plane at 95% C.L. (2 d.o.f, Ax^ = 5.99). In this 
analysis, Stl and St-r are fixed to zero, while the normalization parameters are marginalized. The 
solid red curves indicate the current bounds, while the green, blue, and cyan curves indicate what 
the bounds would be with reduced uncertainty in the ''Be signal normalization (left-panel), and 
reduced ^^Kr background (right-panel). The area outside each contour is excluded. The bound 
from Ref. [16] is plotted in black for comparison. 



If we can reduce this uncertainty from 7% to, say, 3% then for SeL and EeR, the limits can 
be improved roughly by a factor of 1.5. For EtL, we find two disjoint regions if we can go 
below 3%. Such an improvement must first follow the resolution of the solar metalicity 
problem we alluded to in section 4, and further improvements in nuclear cross section 
measurements. 

7.3 Constraints in the {SeL-^eR) and (erL-Srif ) planes 

Let us now turn to 2D constraints in the left-right coupling plane of each flavor. In Fig. 6, 
we show the allowed regions in the (eeL-^ei?) plane at 95% C.L. (2 d.o.f, Ax^ = 5.99), 
obtained with the r NSI parameters set to zero, £rL = ^tR = 0, and the normalization 
parameters marginalized. On both panels, the solid red curves indicate the current bounds 
with 153.6 ton-year fiducial exposure, 7% uncertainty in the normalization of the ''Be signal, 
and with the full *^Kr background. On the left-hand panel, three more curves indicate the 
bounds assuming three different uncertainty levels (see legends in the figure) in the "^Be sig- 
nal normalization. On the right-hand panel, three more curves indicate the bounds with 
three different assumptions (see legends in the figure) on the amount of ®^Kr background. 
For comparison, we plot the combined solar -|-KamLAND bound obtained in Ref. [16] in 
black. In Fig. 7, we show the same for the [EtL-^tr) plane. For comparison, the bound 
based on the LEP 'neutrino counting' data [17] is plotted in black. 

Here, we find that the Borexino bounds are comparable in size to those from Ref. [16] 
and [17], but occupy slightly different regions in the 2D parameter space. Thus, the overlap 
region in 2D is smaller than either Borexino, or the reference bounds alone. While the 
expected improvements to the bounds due to reduced normalization uncertainty and ^^Kr 
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Figure 7. Allowed regions in the ErL-^rR plane at 95% C.L. (2 d.o.f, Ax^ — 5.99). In this 
analysis, e^L and e^R are fixed to zero, while the normalization parameters are marginalized. The 
solid red curves indicate the current bounds, while the green, blue, and cyan curves indicate what 
the bounds would be with reduced uncertainty in the ""Be signal normalization (left-panel), and 
reduced ^^Kr background (right-panel). The area outside each contour is excluded. The bound 
from Ref. [17] is plotted in black for comparison. 

background are modest for the r-parameters, they can be considerable for the electron 
parameters, particularly if the *^Kr background can be reduced to 10%. 

8 Summary 

We have used the 153.6 ton-year fiducial exposure data from Borexino to place bounds on 
the flavor-diagonal NSI parameters SeL, £eR, ^tL, and ErR taking into account the ±7% 
uncertainty in the ''Be solar neutrino flux, $7^^^, and the background from ^^Kr and ^^'^Bi 
/3-decay. The uncertainty in sin^ 623 was not considered. The resulting one NSI parameter 
at a time bounds are listed in Table 2. 2D bounds in the (eeL-^eij) and (erL-eri?) planes 
are shown in figures 6 and 7. They are comparable to existing bounds in Ref. [16] and [17]. 

One caveat to our result is that due to the neglect of the uncertainty in sin^^23, 
the uncertainties on EtL/r are underestimated. A resolution of the disagreement between 
Refs. [53] and [54] is crucial if we are to derive more robust bounds. 

Further improvements in the bounds would require reductions in the *^Kr background, 
which is already underway in Borexino Phase II, and in the uncertainty of the ''Be solar 
neutrino flux. The latter may be achieved by the resolution of the solar metalicity problem, 
and improvements in the relevant cross section measurements by future experiments such 
as DIANA [59]. 
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